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Abstract 

CN I The path integral representation for the Renyi entanglement entropies of integer index 

JT^ I n implies these information measures define operator correlation functions in QFT. We 

^vq ' analyze whether the limit n — )• 1, corresponding to the entanglement entropy, can also be 



O 



^ 



represented in terms of a path integral with insertions on the region's boundary, at first 

order in n — 1. This conjecture has been used in the literature in several occasions, and 

specially in an attempt to prove the Ryu-Takayanagi holographic entanglement entropy 

formula. We show it leads to conditional positivity of the entropy correlation matrices, 

which is equivalent to an infinite series of polynomial inequalities for the entropies in QFT 

or the areas of minimal surfaces representing the entanglement entropy in the AdS-CFT 

context. We check these inequalities in several examples. No counterexample is found in 

o ' the few known exact results for the entanglement entropy in QFT. The inequalities are also 

^ ' remarkable satisfied for several classes of minimal surfaces but we find counterexamples 

corresponding to more complicated geometries. We develop some analytic tools to test 

the inequalities, and as a byproduct, we show that positivity for the correlation functions 

^1 is a local property when supplemented with analyticity. We also review general aspects 

C^ ■ of positivity for large N theories and Wilson loops in AdS-CFT. 

o 
o 

• '• 1 Introduction 

o 

CN ' The positivity of the Hilbert space scalar product gives place to an infinite series of inequalities 
for the correlators of operators in quantum mechanics. Given a family Oi, i = l,...,m, of 
operators, a state |0), and arbitrary complex numbers Aj, we have 

X 

Y, Ko. E ^^^^ io) = E ^*^^- (01^1^.10) > • (1) 

Because the Aj are arbitrary, the matrix of correlators has to be positive definite. This is 
equivalent to 

det({(0|OjO,|0)}„=i...^)>0, (2) 

for any collection of m operators. 

When the starting point in the description of the theory is given by some numerical func- 
tions, these inequalities become specially relevant. In fact, this property allows us to recover 
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the Hilbert space from vacuum expectation values, and plays a central role in Wightman re- 
construction theorem in quantum field theory (QFT) [1]. 

A remarkable case of a function unexpectedly satisfying inequalities analogous to (|2]) is given 
by the traces of integer powers of the local density matrix tr(py), n = 2, 3, .... Here py is the 
vacuum state density matrix |0)(0| reduced to a region V of the space by tracing over the degrees 
of freedom outside V . It is quite surprising these statistical measures give place to operator 
correlators. The operator behind the identification with correlation functions for tr py.y. is 
not {pViVj)^~^, which nevertheless produces the correct expectation values (0 |(p\/iV,)"~^| 0) = 
tipy.y.. Here {pviv)"'~^ is not an operator in the usual sense since it depends on the state. More 

importantly, inequalities such as ([2]) further require that tipyy, = {Oy Oy ), for operators Oy 
and Oy localized in Vi and Vj, in some theory, not necessarily the original one in which p is 

calculated. This is not the case of py~y ^ Pv"~Vy~^- 

An operator expectation value interpretation of tipy follows considering the n-replicated 
model. This construction is valid in any quantum mechanical system. Let us consider a system 
with global state |0), Hilbert space "H = 'Hv<^'H-v, and let pv = tr-^_-^ |0) (0|. In the replicated 
model we take n identical copies of this system. The global state is lO*-"^) = ®"=i |0j), and let a 
basis of the global Hilbert space of the i*^ copy be {|e^) <S) 1/^)}, where {|e^)} is a basis for the 
Hilbert space Tiy and {|/^)} is a basis for the one corresponding to the complementary region 
7i''_y. Then, defining the unitary operator Oy acting on Tiy 

it follows that 

(0(")|C{7^|0(")) =trp^. (4) 

We also have the splitting relation Cy^uvi ~ ^Vi^vJ ^'^^ non-intersecting Vi and V2. A lattice 
version of Oy has been used in [2], where these are called swap operators. Notice Oy does 
not depend on the chosen basis for the Hilbert spaces in ([3]). 

The replica trick was first introduced in this context to give an expression of the trace trpy 
in the path integral formulation. This is given by the partition function in a n-replicated space 
identified along the n copies of V in such a way one goes from a copy to the other crossing 
V. The complete manifold is fiat everywhere except at the boundaries of V where it displays 
a conical singularity of angle 27rn. For a detailed description of this replica trick see [31 SI Ej . 
These partition functions in non trivial manifolds have been interpreted in the two dimensional 
case in terms of field operator expectation values in [5l [6] . These are called twisting operators, 
and are defined to do the job of imposing the correct boundary conditions. 

It is convenient for later use to express these traces in terms of the "entanglement" Renyi 
entropies SniV) = —{n — 1)~^ log(trpy), such that 

tT{pl) = e-("-i)^"(^) . (5) 

Allowing for arbitrary real index n, the limit lim„^i SniV) = S{V) = — trpy logpy, gives the 
entanglement entropy corresponding to V. 

In the Euclidean version of QFT the positivity inequalities analogous to ([2]) are called 
refiection positivity inequalities [7]. Specifically, for the integer n Renyi entropies these later 
write 

det ({trp;:,^^,},,,=i...^) = det ({e^("^i)^"(^'^^)},,,=i...™) > 0. (6) 
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Figure 1: (a) Two regions V and V in Euclidean space. One is the time-reflected image of the 
other, (b) If the regions Vi and Vi are close to the x^ = plane, the common surface near 
this plane should give an overall numerical factor which drops out from the inequalities. This 
configuration gives information about the set formed by the union ^ = ViVi in the limit of no 
separation. 



Here Vi, i = l...m is a collection of co-dimension one sets included in the half-space of positive 
Euclidean time, and Vj is the Euclidean time-reflected region corresponding to Vj. The region 
ViVj is just the union of Vi and Vj (see figure [I]). A generalized version of these inequalities 
valid for any quantum mechanical system has been proved in [8] (in a real time formulation of 
reflection positivity). 

A short and direct proof of the inequalities (^ follows directly in the usual way reflection 
positivity is proved with the path integral. This involves splitting the path integral for positive 
and negative time. Writing 0"*" and 0~ for the restrictions of the field (j) to positive and negative 
Euclidean time, and considering a family of regions Vi in the positive Euclidean time half-space, 
we have 
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These are integrals in n replicated space, and the superscript V means the appropriate boundary 
conditions take place on V. The normalization factor A^ = j Vcjye'^^'^^ is the one copy path 
integral without boundary conditions on V . The fields 0o; 0? 0^ and 0~ have n independent 
components. 0o(^) is the common value of 0"'"(t,x) and 0~(t, x) for t = 0. Provided the 
Euclidean action has the time reflection symmetry S'[0(t,x)] = S'[0(— t,x)]*, a change of field 
variables 0(t, x) — )■ (j){—t,x) in the path integral shows the two terms in the brackets are 
conjugate of each other, and the result is positive. 

In QFT the operators Oy are localized in V, but they do not really depend on the shape 
of V but only on its boundary dV. Different V with the same boundary give place to the same 
boundary conditions of the path integral. In the real time interpretation this is just causality: 
pv and pv are the same if the spatial sets V and V have the same boundary because in that 



case these are equivalent Cauchy surfaces [9]. A unitary transformation between the variables 
on these two surfaces amounts to an irrelevant change of basis in the Renyi operators ([3]). 

The hypersurfaces V have two possible sides. Crossing V from one side, the change of 
copies in the replicated manifold is in the sequence 1 — ?> 2 — ?> 3... or in the opposite direction 
1^2^ 3..., according to the chosen orientation for V. If the orientation is determined by a 
normal vector to the surface V, in the definition of V it is assumed that the orientation vector 
has changed by a time reflection. We can take advantage of this to obtain inequalities for single 
component sets V. Arranging regions Vi close to the surface t = 0, V^Vi approaches a one 
component region V (see figure lb). The sides along the t = plane have opposite orientation 
and their contribution should decouple from the inequality i^ because the nearby operators 
should have an OPE proportional to the identity. 

The entropy is obtained from the Renyi entropies in the limit n = 1. This requires an 
analytic extension in n from the knowledge of Sn for n integer. This extension is unique under 
reasonable assumptions (see [SI [ID])- However, in general there is no path integral representation 
of Sn for general n. The naive approach would be imposing boundary conditions in a manifold, 
with conical angle 27rn, for non integer n, in the boundary of V, and such that for integer n 
one recovers the above mentioned replicated manifolds. Unfortunately this cannot be generally 
achieved using manifolds which are fiat outside dV. A notable exception is the case where 
there is a "rotational" symmetry with the angle 2?™ as the parameter of a one dimensional 
symmetry group (see [TT] for a nice account of these problems). This is the case of V being 
the half space in a general QFT, or a sphere in a conformal QFT. In curved space-times this 
symmetry appears typically when the boundary of l^ is a Killing horizon. 

Another example in which thinking in terms of small deficit angle in the path integral is 
justified is Solodukhin's formula for the logarithmically divergent term in the entropy for con- 
formal theories and general regions in four dimensions [12] . This is given as a linear combination 
of two dimensionless, conformally invariant, and local terms, which are integrals of polynomials 
of the extrinsic curvature on dV. This general geometric structure of the coefficient of the 
logarithmic term can be ascribed on general grounds to the nature of the ultraviolet diver- 
gences in the entropy. The coefficients of these geometrical terms are linear in the two trace 
anomaly coefficients of the CFT. The trace anomaly coefficients are defined for smooth metrics, 
hence it is apparent the necessity of the limit of small conical defects to establish a connection 
with the entanglement entropy. For the case where l^ is a sphere, the coefficient turns out to 
be proportional to the coefficient of the Euler density in the trace anomaly. This has been 
later proved to be the case in any dimensions using purely QFT arguments, which rely on the 
existence of a rotational (conformal) symmetry for the sphere [TH [T3] (see also [S]). For a 
cylinder in flat space the other anomaly coefficient is selected in Solodukhin's formula, though 
there is no additional rotational symmetry. However, the logarithmic terms are both local and 
conformally invariant, and in order to compute them a conformal transformation can be used 
to map the original geometry into another one in curved space, in which there is no extrinsic 
curvature, and there is a rotational symmetry around the boundary (at least locally)lj. The 
result for the cylinder has also been found by holographic calculations of entanglement entropy 
for theories with higher derivative gravity actions [15] , and confirmed in numerical calculations 
of the logarithmic term for free massless scalars and fermions in four dimensions [16] . 

In more general cases, where there is no path integral representation for non integer n, 
one may still think, as a working hypothesis, that for infinitesimal deficit angles 2-iT{n — 1), 
appropriate for the limit corresponding to the entropy, it might be possible to use the path 
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Figure 2: Two regions Vi and V2, their time reflected conjugate ones Vi, V2, and the minimal 
surfaces corresponding to the four combinations ViVj. The surfaces corresponding to V1V2 
and V2V1 cut at a surface f2 at x^ = 0, because the points at x^ = are invariant under 
time reflection. The surface Q divides these two surfaces in four pieces, and application of 
the triangle inequality for the minimal areas leads to f lT2|) . The same argument (though not 
involving Euclidean time-reflection symmetry) was given in [19] to prove strong subadditivity 
of the holographic entanglement entropy for regions lying in a single hyperplane. 



integral formula, and trust it to first order in (n — 1). The conical curvature is then placed by 
insertion of an infinitesimal curvature tensor supported on the boundary of V. This construction 
has been used in the literature in several occasions, though up to the present there is no clear 
deductive way to this representation for the entropy in a general region. 

A notable case where the representation of the entanglement entropy in terms of a path 
integral in a manifold with infinitesimal conical singularities is used is Fursaev's proof of the 
holographic entanglement entropy [17], as given by the Ryu-Takayanagi ansatz [18]. This ansatz 
gives the entanglement entropy in V in the AdS boundary QFT as (4G)~^ times the area of a 
minimal surface in the bulk AdS space, with boundary coinciding with dV, 

S{V) = {4G)-'mme^=evA{^) . (8) 

We often call the area of the minimal surface directly as A{dV). This holographic formula for 
the entropy has passed several tests, including giving the correct coefficient of the logarithmi- 
cally divergent term for the entropy of spheres in any dimension [Ill[l3], the strong subadditive 
inequality for sets in a single spatial hyperplane [19] , and the thermal entropy of the boundary 
in presence of a bulk black hole [20]. However, the missing link between the replicated space 
for integer n and the small angle hypothesis sheds some doubts on one of the assumptions of 
Fursaev's proof for general geometries, as discussed in [TH fll\ . 

Thinking in terms of the positivity inequalities, the idea of the small deficit angle path 
integral representation, with insertions at the boundary of V is far from being innocent, and 
implies a large body of non trivial relations. We have seen above that e~("~^)'^"(^) obeys the 
correlator inequalities due to its path integral representation. The assumption of the small 
deficit angle path integral representation can be then be rephrased (and generalized) as that 
Q-(n-i)Sn(V) Q]3gyg f^i^Q positivity inequalities at first order in n — 1 as n — ?■ 1. This is equivalent 
to the matrices Sij — {n — l)S{ViVj) being positive definite for small enough (n — 1). As explained 
in section 2 this is the mathematical requirement for the matrices —S{ViVj) to be conditionally 
positive definite, or equivalently, the matrices of the exponentials ^-^^i^i^j)^ with any A > 0, to 



be infinitely divisible. Expanding the determinants in (|6]) for small (n — 1) we get 

det{{S{V^j+,) + S{V,^iV,) - S{V^,) - 5(l^,+iF,+i)K,=i„„„_i) > . (9) 

This is an infinite series of polynomial inequalities in the entropy, one for each polynomial 
order. If these conditional positive inequalities where correct the entropy itself would have an 
associated "almost" operator, such that S{V) = {0{V)), ^(^1^2) = {0{Vi)0{V2)). 0{V) is 
an "almost" operator, in the sense that positivity inequalities are not be satisfied by 0(y) but 
they are fulfilled by 0{V) + cy, with cy a numerical function. This function generally has to 
be taken to be non translational invariant. An example of an almost operator in this sense is 
given by the free massless scalar field in two dimensions, as will be review later. 

The entanglement entropies have ultraviolet divergences which might in principle obscure 
the physical content of these inequalities. Correspondingly, the minimal areas in AdS space 
diverge due to the infinite volume of the surfaces as they approach the asymptotic boundary. 
However, these divergences, being local and extensive along the boundaries dV, cancel out in 
the generic matrix element in (Q. The inequalities can be equivalently written in terms of the 
mutual information I{U, V) = S{U) + S(y) — S{UV), which is finite, as 

det{{I{V,, V,) + /(y.+i, V,+i) - I{V,, \/,+i) - /(y.+i, y,)}„=i,...^„i) > . (10) 

The divergences can also be eliminated from ([6]) by computing the well defined connected 
correlators 

^(n-l)I„iV,,Vj) ^ ^{n-l)Snm^in-l)S„{V,)^-{n-l)S„{V,Vj) ^ (^^^-^ 

with In{U, V) = Sn{U) + Sn(y) — Sn{UV). Multiplying the matrix elements in (^ by the factor 
^(n-i)s„(Vi)^(n-i)s„{Vj) ^^Qgg j^q^ change the positive definite character of the matrix. 
The linear inequality in ©(with m = 2) writes 

25(^1^2) > S{ViVi) + S{V2V2) . (12) 

It also holds for the integer index Renyi entropies since it coincides with the m = 2 case in ([6]). 
When Vi C V2 are on the same hyperplane the inequality f lT2|) coincides with some specially 
symmetric case of the strong subadditive inequality of the entropy. This last inequality has 
been discussed in the context of minimal areas [22] • The reason for the validity of flT^ for min- 
imal areas is very simple, and was first discussed in [12] in relation to the strong subadditive 
inequality for coplanar regions. It is due to the triangle inequality satisfied by the minimum of 
the area functional (see figure |2]). However, in the present case, the time-symmetric position 
of the regions in (IT2|) is important for the validity of the inequality, since it enforces the rele- 
vant minimal surfaces cut each other in order to apply the triangle inequality. Further linear 
inequalities for the entropy in AdS-CFT have been discussed in [23]. 

We hope that learning whether the m > 2 non-linear inequalities in ([9]) hold will teach 
us something about the fundamentation of Fursaev's proof and the Ryu-Takayanagi proposal 
for the holographic entanglement entropy. In this paper we take an "experimental" attitude 
and check these inequalities in all (admittedly few) known exact calculations of entanglement 
entropy in QFT. So far, the entanglement entropy can be calculated exactly mostly for free 
fields in simple geometries and these are the cases we have checked. We use either numerical 
evaluation or analytical tools. Up to what we have checked, we could not find a counterexample 
to these inequalities in the QFT examples. 

Then we check them for minimal areas representing holographic entanglement entropies. 
We find the inequalities are remarkably satisfied in several non trivial cases, though we also 



find some counterexamples. We distinguish two types of counterexamples, the ones in which the 
minimal areas in the correlation matrix experiment a phase transition as a function of the shape 
of the region, and the ones in which they do not. These phase transitions are related to the fact 
that the boundary theory is evaluated in a large N and large coupling constant limit, and are 
not expected to occur in the entanglement entropy for ordinary field theories with finite number 
of degrees of freedom. The counterexamples without phase transition violate the "infinitesimal" 
inequalities obtained from (Q for infinitesimally displaced regions. The infinitesimal inequalities 
involve an infinite series of inequalities for the function derivatives of arbitrarily high order. 
Remarkably, we show the positivity relations of correlation functions can be integrated from 
the infinitesimal inequalities to yield the finite inequalities in regions of analyticity for the 
correlators. Hence, when the infinitesimal inequalities hold in the neighborhood of a fixed region 
V, they automatically hold for all the domain of analyticity. These domains of analyticity are 
broken by phase transitions in the entropy. Our results give a negative answer for the general 
validity of conditional positivity for the holographic entanglement entropy. 

The theme of minimal surfaces in AdS establishes a natural connection with Wilson loop 
operators. Maldacena's correspondence gives a geometric prescription for the correlators of 
Wilson loops in a nontrivial interacting theory, the large N limit of A/" = 4 SU{N) super- Yang- 
Mills gauge theories in four dimensions [23]. In Euclidean space and in the the large t'Hooft 
coupling limit A ^ 1 we have 

(iy(C)) = /o(A,C)e--^-^(^), (13) 

where C is the (single component) loop, and A{C) is the minimal area of a bulk two-dimensional 
surface whose boundary coincides with C. This is independent of A^ at leading order, and 
/(A, C) is a non exponential correction. It is then interesting to look at the way the quantum 
mechanical positivity property ([2]) is realized for Wilson loops in this context. 

In the next section we review this issue. We start discussing positivity in the broader 
context of the semi-classical and in large N limits. We recall large N Wilson loops are "classical 
operators" in the interpretation of large N theories as classical limits. We show some exponential 
of classical operators lead to infinite divisible correlation matrices. This is analogous to the 
central limit of probability theory. We then discuss large N Wilson loops highlighting the deep 
differences between the entropy and Wilson loop behaviors as expansions in A^ and A. We show 
the linear triangle inequality for the minimal areas is enough to satisfy positivity at the leading 
order in large A^ and A, and no higher order polynomial inequalities are in principle required 
to hold for the minimal areas for Wilson loops. 

In section 3 we prove positivity properties spread on the domain of analyticity of correlation 
functions. In section 4 we analyze known exact examples of entanglement entropy in QFT and 
check conditional positivity. In section 5 we check the geometrical inequalities (|9]) for exact 
solutions of minimal surfaces in AdS. Finally, we end with a discussion of the results. 

2 Positivity in large N and semi-classical limits 

2.1 Positivity in the semi-classical limit 

The large N limits can be understood in the same terms as the classical limit of quantum 
mechanics [25]. Then we first look at infinite divisibility in the semi classical limit. 



Consider evaluating the expectation value in ordinary quantum mechanics in the Euclidean 
framework 

{q{-ti)q{t2)) = , ^ , . s[.(.)i • (14) 

J Vq{t) e h 

By reflection positivity this has to be a positive definite kernel for positive ti,t2 > 0. Let the 
minimum of the action be attained for some value q{t) = qo- In the /i — ?■ limit the action can 
be expanded around the classical minimum. Writing the deviations from the classical value as 

^i't) = q(t) ~ qo "we have 

S = So + - dhdt2x{ti)K{ti,t2)x{t2) + 0{x{tf), (15) 






Sq{ti)6q{t2 



(16) 

Q{t)=qo 



The most relevant part of the functional integral in the limit of small h comes from values of 
the fluctuations x{t) ~ \/h. Then we have 

{q{-ti)q{t2)) = ql + hK~\-t^M) + O{0^) . (17) 

Positivity for this correlator is trivial to the zero order, but to first order it implies conditional 
positivity of the kernel ^"^(—^1,^2) 

dtidt2a{tiya{t2)K^\-ti,t2)>0, dta{t) = 0. (18) 

The condition on the possible test functions in (llSp is necessary in order to have an inequality 
which does not involve the constant zero order term in the correlator flT7|) . 
In more generality, a matrix bij is said to be conditionally positive if 

x*bijXj>0, "^Xi = 0. (19) 

i 

That is, it is positive for any vector Xi satisfying Xli^* ~ 0- ^^ ^^^ language of QFT b^j is 
almost a correlator, since it satisfies the positivity inequalities except for a one dimensional 
subspace. Equivalently, b^j is conditionally positive if there is a vector ^j such that 

k,+^, + ^, (20) 

is positive definite. Another equivalent condition is that the matrix [22] 

bij + bi+ij+i - bi+ij - bij+i i,j = l...m-l (21) 

is positive definite. 

Another aspect of positivity in the semi-classical limit is shown using correlation functions 
of delta function operators 6{q{ti) — qi). These force the paths in the functional integral to 
pass through the point {ti,qi). As a result, in the classical limit, the correlation functions are 
dominated by exponentials of the classical action 

(SiqiU) - qi)S{q{t,) - q,)) = ce-^-«*-''')'(*-''^))/''(l + Oik)) . (22) 



Then the reflection positivity inequahty for two points boils down to the triangle inequality 
for the classical action functional, in a way analogous to the one shown in figure ([2]) for the 
minimal surfaces. We have 



e 



-S'ci({-ti,gi),(ti,gi))/ri. p~SM-ti,qi),it2,q2))/fi 



^^M e-'5'cl((-«2.92).(tl.9l))M (,-SM-t2,q2),it2,q2))/h 1>0, (23) 

which is equivalent to 

5ci((-ti, gi), (ti, gi)) + 5ci((-t2, ga), (^2, ^2)) < 2S,i{{-h, gi), (ta, g2)) • (24) 

The triangle inequality is immediate for a real Euclidean action, containing at most first deriva- 
tives of the fields. There is a possible failure of the triangle inequalities for higher derivative 
theories due to boundary terms in non smooth surfaces which are required by the proof. This 
is related to the failure of these theories to be unitary. 

It is not difficult to convince oneself that the triangle inequality for the euclidean action is 
all what is required for the positivity of correlator matrices involving several delta function op- 
erators, to leading exponential order in h. That is, in the fi — )■ limit, the inequalities for higher 
order determinants for delta function operators are satisfied provided the two-dimensional de- 
terminant inequality holds, since it is the only one homogeneous in h. 

2.2 Infinite divisibility, the central limit and free fields 

In this subsection we comment on the relation between conditional positivity and another 
property called infinite divisibility, which is related to the central limit theorem. 

Infinite divisible property is inspired by Schur's theorem for positive definite matrices: If 
{aij} is positive definite (it is hermitian and all its eigenvalues are positive) then the matrix of 
the integer powers of the elements {(o-jj)^} is also positive definite. A positive definite matrix 
{ttij} = {cfj} which is equal to the entry- wise A^-power of another positive definite matrix 
{cij} for any A^ (here {cij} depends on A^) is called infinite divisible [27]. In QFT, the infinite 
divisibility of a correlator, means that its fractional powers also satisfy the positivity conditions 
a correlator ought to satisfy. Then, these fractional powers can be thought as correlators for 
operators in another theory, non necessarily the original one. 

It is not difficult to show that the positive definite matrix {aij}, i,j = l...m, is infinite 
divisible if and only if the matrix of the logarithms of its entries 

bij = log(aij) (25) 

is conditional positive |26] . 

An example of correlators which are powers of other correlators is given by a QFT which is 
the tensor product of identical independent sectors named by 1, ..., A^. Choosing the correlators 
of product operators 

(0|OP...Of )t X Of\..Of^\0) = {0\Ol''>^Of\0)...{0\Of^^Of^\0) . (26) 

Most correlator matrices in QFT are not infinite divisible. In principle, this property should 
be valid only for special type of theories and special types of operators in these theories. Eq. 
(126|) suggests that in order to have infinite divisibility one has to consider a theory in the limit 
of a large number of degrees of freedom, and some special kind of operators. We show in the 
next section this is the case of some exponential operators in large N theories. 

9 



Infinite divisibility is also related, and for the same reason which operates in large N theories, 
to the central limit of probability theory [28]. Take the probability distribution for a variable 
y = -4^{xi + X2 + ■■■ + X]\f) which is proportional to the sum of A^ independent variables Xj, 
with the same distribution P[x). Then the probability distribution Q{y) for y is the A^-fold 
convolution 

Q{y)= dxi...dxN6(y-^{xi + ... + XN)jP{xi)...P{xN). (27) 



The expectation value for the exponential Gqlk) = (e*'^'^) (called the characteristic function of 
the distribution) is 



Ggik) = j dye^'yQiy) = {e'vWy = Gp{k/VNy . (28) 

Then, the infinite divisible case for the characteristic functions corresponds exactly to proba- 
bility distributions which arise from an infinite sum of variables. Under the condition that the 
variance of x is finite, the limit of the sum of infinitely many variables can only be the Gaussian 
distribution, which of course has infinite divisible characteristic function. This is the central 
limit theorem. Different central limit distributions are also possible when the variance is non 
finite. 

Notice that the characteristic functions satisfy positivity relations which are similar to the 
ones corresponding to two point correlators in QFT: 

X*GQ{h - kj)X, = I dxQ{x){\e-'^y{\je-'VN-) > 0. (29) 

Thus, Gglki — kj) are positive definite matrices, and f l28|) means these write as entry-wise A^ 
powers of other positive definite matrices. 

Gaussian distributions are analogous to free fields, and the proof of infinite divisibility for 
the expectation value of exponential operators is direct. We have for a free field (f){x) 

/Q|g«/'^^fc(a:)0(a:)|g\ _ ^-^ J dxdy k{x)g(x-y)k[y) /ggN 

where g{x — y) = (O|0(x)0(i/)|O) is the two point function, and k{x) is an arbitrary function. 
Then 

^Q|gi/dxfc(x)<)!.(x)|Q^ = (0|e*-^'^^ vi'^^^^lO)^ = (o|e*-^'^^'"^''^^^'^'^^^"^-+'^'^^^^^|0) . (31) 

The first equation shows the expectation value of an exponential operator as a power of the 
expectation value of a different exponential operator in the same theory, and the second equality 
shows it as an expectation value in a different theory, which is one of A^ independent free fields 
0i(x),...,0Ar(x), with the same two point function. 

Any exponential of a free field gives an operator with infinite divisible correlators. Another 
example is the Wilson loop operator e'^^3"i^''^t^ for the free electromagnetic field. Notice that 
the two point function of the free field g{x — y) has to be only conditionally positive, and 
this allows for operators which are not strictly speaking quantum fields. This is the case of a 
massless scalar field in two dimensions, whose two point function {0\(j){x)(j){y)\0) ~ — log |x — y| 
cannot be positive for all range of coordinates x, y. In section 5 below we show — log \x — y\ is 
a conditionally positive correlation function. 
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Examples of infinite divisible operators are given by the exponentials e v^ of classical 
operators divided by \/h is the classical limit. We have 



F[q{t)] _ F[qo] , f ^^^ 6F[q{t)] 



Vh Vh J Sq{t 



<h) ,o('^]. (32) 



<?(t)='7o 



Vh \ Vh 



Higher order deviations from the classical value in F[q{t)] or the action are suppressed for small 
h. In this limit, the operator will behave as an exponential of an operator linear in a Gaussian 
field, and we expect infinite divisibility. We have 



(e^) = e4re^/'^*^^*^'^(*^)(''"^)(*^'*^)'^(*^) (l + 0(v^)) , (33) 



where J(t) = -j-J^ . Given a collection of functionals Fi[q{t)] with support on the positive 
time t > region, this later requires the positivity of the matrix (to leading order in h) 

Fjlqi't)]" Fj[q{t)] Fjlgp]* Fjlgp] 



e 



^Tir^e^T^) = e^^e^^e^^'^*^'^*^(-^>*(-*^)+-^^(*i»('^"')(*i'*^)('^^(-*^)+-^^(*^^) . (34) 



The positivity of the matrix in ( 15^ is equivalent to the conditional positivity of the kernel 
(K~^)(ti, ^2)- We have seen this also follows from the positivity of the two point function to 
first order in h (TTSj) . 

Using functional integrals in phase space, one can also show in the same fashion the infinite 

Flq,p] 

divisibility of operators of the form e v^ , provided the functional F[q^p] is well behaved in the 
h -^ limit, and that F[q,p\ is such that the functional integral converges. The operators such 
as F[q,p] which in the limit 7i — ?■ can be expressed as a function in phase space (technically, 
they have a well defined limit of expectation values in any coherent state) are called classical 

Flq,p] 

operators in [25]. Thus, the infinite divisible operators e ^ are not classical operators. They 
are half way between classical operators of the form e^^'''^'^ and coherent operators of the form 
e~i^ . The vacuum expectation values of operators of the form e~ft^ produce much greater 
excitations in the /i — )■ limit, and test terms in the action which are of higher order that the 
Gaussian term. Thus, in general there is no infinite divisibility for these operators. 

In the large N limit of interacting theories [29] a result on infinite divisibility analogous to 
the one for free fields depends on choosing adequate operators that test only Gaussian fields 
fluctuations in leading approximation. 

2.3 Large N theories as classical limits 

According to the general scheme of [23] the large N limits can be thought as classical limits. 
This is done by adequately choosing a set of classical variables in the large N theory, and 
identifying the small parameter x, which plays a role analogous to h in the classical limit, 
with some inverse power of N. The expectation values are then computed with a path integral 

over the classical variables, with a weight given by e ^, S^ being a function of the classical 
variables, regular in the x ~^ limit. 

Then, as in the classical limit discussed previously, positivity implies for example the condi- 
tional positivity of the connected correlators of classical operators. Equivalently, it also implies 
the infinite divisibility of operators which are the exponentials of the classical operators divided 
by y/x- While the computation of the action in terms of the classical variables, or evaluating 
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its minimum, can be very difficult in specific problems, we do not need this information here. 
We only need to identify the classical operators. This has been done already in [25] for a large 
class of models. 

The general scheme (for more details see [23]) starts by defining a coherent group G and a 
unitary representation G^ for each x- The group G depends on the problem. Choosing a base 
state |0) , the coherent states are defined as the states \u) = u |0) for any u G G^. These form 
an over-complete basis of the Hilbert space. The importance of the coherent states is that their 
overlap in the x ~^ limit is exponentially small, and they can be used to write down a path 
integral. The classical operators can be defined as the ones which have a well defined limit of 

iU \A\ U')y , . 

lim V }'"" . 35 

The interest to us here is that a generating set of the classical operators is obtained from the Lie 
algebra of the coherence group. If A belongs to the Lie algebra of G^, then (x^.) is a classical 
operator. The classical operators are generated by these ones. From this, natural candidates 
to infinite divisible operators are of the form e*^^^. 

The example of the classical limit can be useful here. The coherence group is given by 

operators of the form 

^ = ei(ap+feq+c)^ (36) 

The operators generated by the Lie algebra elements times h are just the algebra of polynomials 
in p, q, which are classical operators. 

Then, in order to find out which operators are infinite divisible in each model we only need 
to know the lie algebra of the coherence group. 

For large N U{N)-gauge theories, which we are interested here consider, following [2S|, a 
lattice model with link variables V^ which are N x N unitary matrices, and a labels a oriented 
link on the lattice. The conjugate momentum E"' is a hermitian matrix, and the commutation 
relations write 






rpa rpP 



^^""4,^.r , (37) 

^5"^(4,E.° - 5uEt,) . (38) 



The Kogut-Susskind Hamiltonian is invariant under local U{N) transformations, and is given 
by 



[ A 5^(^")' - \-' 5^(V^^^ + V^P) j 

\ a p / 



H = N'tT X > {E'^y - X-' > {V'P + V'P) . (39) 



Here the t'Hooft coupling constant A = g'^N is chosen fixed such that the theory has a mean- 
ingful large N limit. The normalized trace is tr = (A^~^)tr, p and p represent a lattice plaquette 
with the two possible orientations. More generally write V^ for the ordered product of the link 
variables along the lattice path F. 

In this case the Lie algebra of the coherence group is infinite dimensional, and the elements 
depend on arbitrary closed loops F on the lattice. It is generated by the operators 

I N^ t~r 1/^ , iN^ (t~r E'^V^ + V^E'') . (40) 

Notice that the Wilson loop operators are precisely of the form trV^ , which are traces of ordered 
products of the link variables along F. The small parameter is here x = ^~^ [2S]- Thus, the 
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Wilson loop is a classical operator. This is the reason it has a N independent expectation value 
(14T]) in the large A^ limit. Thus we expect the connected correlators to be conditionally positive 
in the large A^ limit. 

2.4 Positivity for Wilson loops in AdS-CFT 

The AdS-CFT prescription for the correlation function of Wilson loops is given byn 

{W{C,)W{C2)) = WCi))(l^(C2)) + ^/2(A,Ci,C2)e-^-^(^i^^) + (9(AVA^^), (41) 

This expansion is organized in powers the string coupling X/N ~ Qs and it is assumed the limit 
of large A^ is taken first than the limit of large A. The A{CiC2) is the minimal area of a surface 
connecting Ci and C2, which can be either connected or disconnected depending on the relative 
position of Ci and C2. 

Therefore, as expected from the general properties of the large A^ limit, the connected 
correlators vanish to leading order in A^. This is very different to what is expected for the 
entropies of two separated regions. The first term in the large N expansion of the entropy 
S{ViV2) is expected to be of order A^Mn d = 4 [E] (that is ~ {G^^)R\ with R the AdS 
radius), independent of A (in the large A limit), and proportional to the global minimal area 
for all surfaces with boundary in dVi U dV2, independently of its topology. Otherwise the 
mutual information would be always subleading in A^ with respect to the entropy, and this is 
not possible since when the boundaries of Vi and V2 are near the mutual information must tend 
to twice the entropy. 

Thus, while there is direct competition between surfaces of different genus for the leading 
term in the entanglement entropy, this competition occurs at higher order in A^~^ for correlators 
of Wilson loops. This is because the contribution of surfaces of higher genus to the Wilson loop 
have additional powers of the string coupling. 

Wilson loops satisfy a reflection positivity relation analogous to ([6]), 

det {{{W{Q)W{Cj))h,,=i...^) > . (42) 

Here Ci, i = l...m is a collection of loops included in the half-space of positive Euclidean time, 
and Cj is the Euclidean time-reflected loop corresponding to Cj. The vector tangent to the 
loop indicating the circulation direction gets multiplied by — T, where T is the time reflection 
matrix. 

Taking the large A^ limit of the inequality f H2]) . and using the expression ( I^Tj) for the 
correlators, we have the conditional positivity of the matrix 

/2(A, Cj, Cj)^ ^/\(a(c,)+a(c,)-a(c,c,)) \ /.ON 

/o(A,a)/o(A,C',) /• ^ ' 

In the present case, the stronger condition of positivity holds for this matrix. To leading 
exponential order in A this is equivalent to the positivity of < g-V^AiCiCj) I ^^ ^^ ^j^g infinite 
series of inequalities 

det({e"^-^(^»'^^)|) >0. (44) 

V L J / i,j=l,...,m 



^Wc thank ,Juan Martin Maldacena for clarificatfons on the subject of this section. 
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It is not difficult to see that due to the large A limit all these inequalities collapse to the 
linear inequality for two regions, which is the case m = 2 in dS]), and is the only inequality 
homogeneous in A, 

2A{C^C2) > AiCiCi) + AiC2C2) . (45) 

This is analogous to flT^ . The reason for the validity of fHSl) is again the triangle inequality 
of the minimal areas, the same reason as in the case of holographic entanglement entropy [TH] . 
Therefore, at least in the large A limit, the positivity of Wilson loops holds provided they are 
given by exponentials of minimal areas with large coefficients, without further requirements on 
the nature of the surfaces or the metric. In particular no higher order inequalities such as ([9]) 
are involved. 

According to (H2]) for m = 2, a. general version of (HSj) . the linear reflection positivity 
inequality for the logarithm of the expectation value of the Wilson loops, 

log{W{Q)WiC,)) + log{WiC,)W{C,)) > 2 log{W{a)WiC,)) (46) 

is valid in any gauge theory independently of the large N, large A limit. This was discussed for 
example in [30] . It implies the behavior of the Wilson loop is bounded between a perimeter and 
an area law. A generalization inspired in strong subadditivity has been discussed in [21] for the 
case of Wilson loops which live in the same two-dimensional plane in the AdS-CFT context. 

Since we are considering the large A^ and large A limit we can think in rearranging the 
expansion (14T|) to allow for the large A limit to be taken ffist than the large A^ limit. In this 
case the correlation function of Wilson loops is directly dominated by the exponential of the 
minimal surface, disregarding of the topology of the surface. It is interesting to note that in this 
case we are in presence of a different semi-classical limit, with a different parameter /i — > A~^ 
rather than h — )■ N""^. In this limit the Wilson loop should be considered as a delta function 
operator in the space of loops rather than a classical operator (in analogy with the discussion 
at the end of section 2.1), and the correlators are dominated by the exponential of the classical 
action on the trajectory joining the different points. The inequality fHSjl is then analogous 
to ( 12^ . which follows from the triangle inequality for the classical action. This explains why 
positivity of flH|) holds in addition to conditional positivity of 



3 Positivity and analyticity 

In this section we analyze the interplay between analyticity and positivity properties. We have 
in mind the inequalities for the entropies, but the arguments are generally valid for correlators 
in any QFT. The positivity inequalities in QFT have different formulations. The inequalities 
coming directly from positivity of the scalar product in real time ([2]) always involve the product 
of operators at coinciding points, and the distributional character of the correlation functions 
appears in a fundamental way. We are using here the reflection positivity inequalities in Eu- 
clidean time (where, on the other hand, the minimal area prescription holds in the AdS-CFT 
formulation) which does not involve correlators for coinciding pointc|. 

The conditional positivity of the entropies is given by inequalities ([9]). According to the 
general results of section 2, all the matricial inequalities ([9]) can be summarized as a single 



^There is also a real time formulation of the reflection positivity inequalities which does not involve products 
of operators at coinciding points. These "wedge reflection positivity" inequalities [3 [32] write in the present 
context as ([H]), but the regions Vi are spatial regions included inside the right wedge x^ > \x^\ in Minkowski 
space, and the reflection operation is the wedge reflection, x°' = —x^, x^' = ~x^ , x^' = x^ , j ~ 2...d. 
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relation expressing the conditional positivity of the entropy (minimal area) kernel: 

- f VVi VV2 F[Vi] S{ViV2) F[V2]* > , I VV F[V] = , (47) 

for any functional on regions F[V] with support on the positive time half-space. The measure 
in region space is not relevant in (147|) since it can be absorbed in F. Eq. (Q follows from this 
one by considering F[V] as a sum of localized delta functionals in the different Vi. Eq. (147|) 
follows from (Q by taking a discrete approximation to F[V] as a sum of delta functions. 

In the following, we consider a set of regions described by a finite number of parameters 
which live in an open set D E R'. Let us start with the single parameter case q = 1- Eq. Wf\ 
then writes 

- dxdy(f){x)S{x,y)(t){y)*>0, dx (j){x) = , (48) 



where 0(x) is any test function with support in D, and we have written S{x,y) = S{ViV2), 
for the regions Vi and V2 described by one parameter, x and y respectively (notice y is the 
parameter corresponding to V2, not to V2). In order to fulfill the second condition in fHHj) we 
take (f){x) = dx^{x). Then, integrating by parts we have 

dxdy^{x)h{x,y)^{yy>^, (49) 

where /i(x, y) = —dxdyS{x, y). Therefore, the conditional positivity of S becomes the positivity 
of/i. 

Now, we want to see the implications of (HH]) for coordinates in a very small region in the 
neighborhood of a point x. We assume fi{x,y) is real analytic around a point {x,x). In this 
neighborhood we write any coordinate a.s y = x + e. We use again the discrete version of 
(H9|) involving only a finite number M of points, which follows taking c/? as a sum of localized 
delta-like wave packets. We have 



a, 13=1 a,/3=l m,n=0 

for any M, Aq and e^ with a = 1, ...,M. Given a fixed M, we choose e^ = ee^, with all the 
ia different. We will later take the limit e — > 0. Then, the matrix (eQ,)"/n! for a = 1...M, 
and n = 0...M — 1, is invertible. This follows from the Vandermonde determinant. Hence, for 
a fixed e and e^ we can choose the A^ such that Vm = J2a=i '^a(fa)™'/''^' ^^^ "^ = 0...M — 1 
is any chosen eigenvector of {fi^'^{x,x)}^^~^Q, with unit norm. Lets call (3 the corresponding 

eigenvalue. Since the Aq are at most of order e~^^^~^\ the sum Um = J2a=i ■^ol^a)™'/''^' ^^^ ^^J 
m > M is then of order e at least. We then have 

E ^"^^ E /r"(^'^) ^^ = /5 + E E /r(^^^)^'"^< + 

a,l3=l m,n=Q ' ' m=0 n=M 

00 M—1 00 00 

J2 E /r"(^' ^)^-< + E E /r"(^, ^)^r.< = /? + o(6) > o . (si) 

m=M n=0 m=M n=M 

Taking the limit e — ?> 0, this gives /3 > 0. This is equivalent to say that the matrix /™'"(x,x), 
m, 77, = 0, ..., M — 1 is positive definite. Therefore, we obtain the infinitesimal inequalities 

det({r'"(a;,a;)}^^;io)>0 (52) 
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for all M. Hence, the kernel f^'^{x,x) is positive definite in the space of all m,n. 

Now, for an analytic function, these infinitesimal inequalities imply the ones for finite values 
of the displacements from the point x. Writing (1501) as 



m,n=0 \ a / \^=1 / 

we see this inequality holds automatically once the infinitesimal inequalities are valid. 

The equivalence between the infinitesimal and finite inequalities in the convergence neighbor- 
hood of the Taylor series have a very important consequence. Once the infinitesimal inequalities 
hold at a point x they automatically hold at any other point in a neighborhood, and then the 
infinitesimal inequalities hold in any point inside the convergence radius around x. Regions 
of validity of the infinitesimal inequalities can then be extended in the same way analyticity 
domains are extended. The infinitesimal inequalities have the same domain as the analyticity 
domain. The finite inequalities hold (at least) provided the points in the correlator matrix are 
contained in the convergence radius of the Taylor expansion around some point in this domain. 
Since correlation functions are usually analytic in QFT, the positivity turns out to be quite a 
local property. 

From (152|) we see the infinitesimal inequalities contain derivatives of the function /i of all 
orders. Specifying for the case of a translational invariant entropy which only depends on the 
Euclidean time variable x, S{x,y) = S{x + y), the inequalities write 

/m(x) = det({-^™+"+^(x)}^^;io) > . (54) 

More explicitly the first three inequalities read 

/i(x) = -S"ix)>0, (55) 

Mx) = S{x)"S{x)""-{S{xy"f>0, (56) 

f.,(x) = iS^^\x)f + S^''\x)iS^'\x)f + iS^^\x))^S^^\x) (57) 

-S^^\x){2S^^\x)S^''\x) + S^^\x)S^''\x)) > . 

Writing /o(x) = 1 we have a simple recurrence rule for the successive inequalities, 

/„+i(x)/„_i(x) = /:/„(a;) - (/;(x))^ (58) 

fn{x)>0, n = 0,l,2,... 

All these inequalities are clearly independent. It is not difficult to find solutions of /i(x) > 
which are not solutions of f2{x) > (i.e. S{x) = tanh(x)). 

Perhaps it is interesting to see in a simple example how subtle the inequalities (l58|l can turn 
for /„ with large n. Take for example the function 1/x. Being a correlation function in some 
QFT, gives place to positive definite correlation matrices. Then taking fi{x) = 1/x, x > 0, all 
the inequalities are satisfied. Now consider writing fi{x) = 1/x — a for some positive a. This 
is not any more a positive definite function since for large enough x it changes sign. However, 
for small x the positive definite correlator dominates and one can expect the inequalities are 
satisfied. A simple calculation shows that fn{x) changes sign at x = — . Thus, for any x there 
is a negative fn{x) for large enough n. This is in accordance with the above result on the 
local character of the inequalities for analytic functions: The non positivity of 1/x — a can be 
detected with the infinitesimal inequalities at any x, no matter how small. This also shows 
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that the /„ for large n are sensitive to smaller disturbances. On the other hand, if we consider 
/i(x) = 1/x + a with positive a, this is a positive definite function. It then follows that the 
short distance leading term /i(x) ~ 1/x has to be a positive definite function by itself. This 
shows that in checking the inequalities, it is also useful to check positivity for the leading short 
or long distance approximations. 

When the region depends on more than one variable a similar analysis can be done. Let 
us write collectively these variables as x = {x'^),fi = 1, ...,q. In a similar way as for the one 
dimensional case, the conditional positivity of —S{x, y) implies the positivity of the vectorial 
kernel given by the derivatives of the entropy, fas{x,y) = —d^d^S{x,y). Expanding this in 
power series around a point x we get 



^ 'rlyi (fq^fJ-q 



a,P=l a,/3=l (M)(.) ^'- ^^' ^'- ''''■ 

(59) 
where we have used the multi-index notation (/i) = (/xi, /i2.../Xg), with yUj = 0, 1, 2, .., for i = l...q. 
The strategy is the same as for the one dimensional case. First, we choose the maximum 
derivative order p = max(/ii + /i2 + ... + fiq) which will appear in the infinitesimal inequality. 
This maximum number of derivatives includes a number of different derivative types (/i) given 

(p -\- a \ I P ~^ Q \ 

1331 . Then, we choose the case of M = , and the X"^ 

q J \ q J 

such that 

fT(/x) _ 

is a normalized eigenvector of /^^^ (a^,a^) (in this matrix the derivative types (/i) and {u) 
belong to the first M types). This can be done because it is always possible to choose e\ such 
that the M x M multidimensional Vandermonde matrix (e^)^i...(e|^)^« is non singular |34] . 
Then, the same reasoning as in l^5T\} leads us to conclude that all the eigenvalues of fls {x, x) 
(a qM X qM matrix) are positive, and this matrix is positive definite. These are the differential 
inequalities in the multidimensional case. Again, the differential inequalities imply the finite 
inequalities in the domain of convergence of the Taylor series. This in turn leads to an automatic 
extension of the validity of the differential inequalities to all the domain of analyticity. 

In order to clarify the discussion, let us be more explicit and look at the example of two 

p + q 



v^-' = Y^»-f^--f^ (60) 



variables, q = 2. For each p = 0, 1, 2... we have that a square matrix of order q 

{p + 2){p+ 1) in the derivatives of S{x, y) is positive definite. However, there may be repeated 
lines and columns in fls i^^^) because of the derivatives in a,S. Since a sub-matrix of a 
positive definite matrix is also positive definite, we can eliminate these redundant lines. The 
first case p = involves only the type (yu) = (0,0), and contains, due to the a,S derivatives, 
only first derivatives on each of the entries of S{x, y), 

5(i'0)'(i'0)(x,x) 5(i'0)'(0'i)(x,x) \ 

5(0'i)'(i'0)(x,x) 5(0'i)'(0'i)(x,x) J ■ ^^^ 

The positivity of this matrix is a consequence of the linear inequality 0121) . and always holds 
for a minimal surface which is the true minimum, not only an extreme of the area functional. 
The second infinitesimal inequality p = 1 involves the types (/i) = (0,0), {^) = (1,0) and 
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(/i) = (0, 1), and contains one or two derivatives in each entry 






5(0'l)'(0-2)(x,x) 



S{m,{o,i)(^x,x) 5(i'0)'(^'i)(x,x) 5(i'0)'(2.o)(x,x) 
Sio,imi)(^x,x) ^(0'i)'(i'i)(x,x) ^(0'i)'(2.o)(x,x) 

^(1'1)'(0'1)(X,X) 5(1'1)'(1'1)(X,X) ^(l'l)'(2-0)(a,^2;) 
^(2'0)'(0'1)(X,X) ^(2,0),(1,1)(3,^3,) ^(2,0),(2,0)^3.^3,) 

^(0'2)'(0'l)(x,x) 5(0'2)'(1'1)(X,X) 5(0'2)'(2'0)(x,x) 

(62) 
We recognize the upper 2x2 sub-matrix is just the matrix flM]) . and the positive definiteness 
of dnSD entails the one of flM]) . However, in the positivity of flH2]) . three more inequahties are 
added. Hence, it is not enough to check the positivity of the determinant but to check the one 
of all the eigenvalues is necessary. 

More generally, the analysis of this section can be generalized to the case of infinitesimal 
inequalities for correlator matrices (or conditional positive entropy matrices as we have focused 
the discussion here), of different operators Oi{xi) evaluated in the neighborhood of different 
points Xi (there may be also different number of parameters determining each Xi). The infinites- 
imal inequalities imply the finite ones inside the radius of convergence of the Taylor expansions. 
This provides a further extension of the inequalities, which propagates inside the multidimen- 
sional analyticity domain of the correlator matrices. Basically, once positivity holds at some 
neighborhood of a multi-point Xj, the only obstacles for further expansion have to be sourced 
by singularities which break up analyticity. 

4 Conditional positivity and entanglement entropy in 
QFT 

In this Section we check conditional positivity of (minus) the entanglement entropy in QFT 
using numerical and analytical methods. The available exact results for the entropy in QFT 
are not many, and involve mainly free fields in low dimensions and simple geometries [35] . The 
numerical tests require high precision, typically around 15 — 20 good digits for the /s defined in 
the previous Section. This leaves us with a rather restricted scenario, both from the analytical 
and numerical point of view, and precludes the use of approximations such as the ones obtained 
by putting models in a lattice. 

The studied examples are consistent with conditional positivity of —S{y). However, by the 
above mentioned reasons we are not able to distinguish if this is due to the simple geometries 
used, or the free character of the models. In this sense, we have to mention two facts. The first 
one is that the counterexamples found for infinitesimal inequalities on the minimal surfaces 
in the next Section occur for more complicated geometries for which there are no pure QFT 
calculations of the entanglement entropy. These involve the entropies of two disjoint circular 
regions with different radius. The one parameter functions which are our best QFT cases 
also seem to be conditional positive in the minimal surface case (with the difference that 
only conformal examples are given by the minimal surfaces, while here we have also massive 
examples) . 

The second commentary is that the case of free fields is up to a certain extent special since 
in this case the Renyi entropies are also found to satisfy the conditional positivity inequalities 
while this does not hold for some interacting cases we review below. There is also an indication 
that the traces trp^ with non-integer index n might be positive. If this "fractionality" property 



is correct, extending the positivity to non-integer n, it would naturally lead to infinite divisibility 
of the exponentials of the entropy in free fields. 

However, conditional positivity in the free case does not hold for discrete systems. There is 
a natural generalization of the reflection positive inequalities applicable to discrete systems [S]. 
We have applied it for free discrete fermions (though we are not presenting details of this cal- 
culation in this paper), and have found counterexamples to the inequalities for some numerical 
examples using random density matrices in low dimensional Hilbert spaces (though conditional 
positivity holds for a large fraction of the random examples, specially as the dimensions are 
increased). Hence, if conditional positivity holds for free QFT, the continuum limit has to be 
an important ingredient. 

4.1 One interval in two dimensions: conformal and large mass limits 

The entanglement entropy for an interval of length L in two dimensions and any CFT is given 
by Hi 

S{L) = I log(L/e) , (63) 

where c > is the Virasoro central charge, and e a short distance cutoff. The logarithm is 
proportional to the correlation function of massless scalar fields (0(L)0(O)) ~ — log(L), hence 
it is a conditionally positive function. We give a direct proof of this in the next Section. 

On the opposite extreme to the conformal case, the long distance, large mass leading term 
for the entropy of an interval in a massive (interacting) theory is also known to have a general 
expression given by 

S{L) = -- V Ko{2m,L) + const , (64) 

o ^-^ 

where the sum is over the spectrum of massive particles of the theory with masses rrii, and 
Kq[x) is the Bessel function. This formula has been shown to hold for generic integrable 
models [HIESI, and argued to be valid for any massive two dimensional QFT [57]. For free fields 
a direct calculation is in [38| [39] . Again, we find this entropy is proportional to a correlator. 
This is the one of massive free scalar fields in two dimensions. Hence, it is conditionally positive. 
In fact, this contribution is positive (rather than conditionally positive) since the massive scalar 
is a well defined field operator. 

4.2 One interval in two dimensions: massive free fields 

Still in (1 + 1) dimensions, we now study the case of massive scalar (S) and Dirac (D) free fields 
in a one interval set. The entropy function S{y), with y = mL, is known exactly and given by 
(see also [3S]) 

My) = I db—^:j—-^Sl;'{y), (65) 







sinh(7r6)" 



IT 



Ss{y) = I db—-—S^{y), (66) 

/o cosh(7ro)^ 
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Figure 3: The logarithm of the functions /j, i = 1, ..., 8 (from bottom to top) of the parameter 
y = mL for the fermionic entropy S{y), computed with a series expansion around the origin 
up to y^^ and y^^. The bifurcation between these two approximations at the end of the curves 
indicates the hmit of the region where the series approximation can be trusted. The logarithm 
is real showing these functions are positive. 
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(67) 
(68) 



Equation ( 168|) is a non linear ordinary differential equation of the Painleve type. The difference 
between the fermionic and scalar cases is exclusively due to the boundary conditions on this 
equation, 



uoiy) -^ -smh{b7i)Ki2b{y) as y ^ oo , 

71 

2 
us{y) -^ -cosh{b7i)Ki2b{y) as y ^ oo 

TX 



(69) 
(70) 



The functions —Sb{y) are not conditionally positive, since their large distance approximation 
Ki2b{y) does not satisfy the inequalities. However, our numerical tests for the entropies fl^S]) 
and 066|) using random interval sizes are compatible with conditional positivity. It is not easy 
to attain excellent numerical precision for these integrals except for the case of the Dirac field 
where the solution for the differential equation can be obtained by a series expansion around 
the origin in terms of powers of y and log{y) [35]. The figure ([3]) shows the functions /«, 
i = l...,8, calculated from the fermionic entropy in a region around the origin for the series 
expansion up to orders y^^ and y^^. These functions are positive, what is compatible with 
the conditional positivity of —S. For the scalar case, the expansion for small y starts as 
S{y) ~ I log(y) + I log(log(y)), which also satisfies the inequalities. 

The Renyi entropies are obtained from an analytic extension of Sb for imaginary values of 
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Figure 4: The logarithm of the functions /i, /2, /s, /4, /s (from smaller to higher slope) of the 
parameter y = mL for the fermionic function Sia in the case a = 1/3. The logarithm is real 
showing these functions are positive. Similar results hold for different a and for the scalar field. 
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fe=-(n-2)/2 
(n-l)/2 



(71) 
(72) 



fc=-(n-l)/2 



We checked numerically the differential inequalities fISS]) for the analytically continued functions 
5*4(1, both for scalar and Dirac fields, and for several values of a G (—1/2, 1/2). In all cases we 
find these inequalities are satisfied. This points to conditionally positivity of these functions 
and of the Renyi entropies —Sn- The figure (j4]) shows a particular case a = 1/3. As shown 
in the figure, all the curves fi look very similar. This may indicate the conditional positivity 
in this case could be a consequence of the interplay between the recurrence rules (1581) for the 
inequalities and some recurrence relation on the Painleve differential equation (l^Sj) . 

The same rules fl65|) . f l66|) . ( FTTj) and fl72|) relates the entropy and Renyi entropies for free 
fields of any region and space-time dimensions, perhaps pointing to the generality of this rela- 
tion between conditional positivity for the entropies, Renyi entropies, and fractionality in free 
theories. 

It is interesting to note that the exponentials of the Painleve-related functions S^a and the 
entropies, would then belong to the class of positive definite infinite divisible functions. Infinite 
divisible functions which are also probability distributions have been classified in relation to 
the central limit theorem in probability theory [28] . 

4.3 Multiple intervals for a free massless fermion in two dimensions 

Our next example is given by the massless fermion in two dimensions [38]. This is, so far, the 
only known complete entropy function. The entanglement entropy corresponding to a set V 
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formed by p disjoint intervals (a.j, 6j) on a spatial line, with ai < bi < flj+i, is 

S{{ai,bi)...{ap,bp)) = - f ^\og\ai - bj\ - ^log ja^ - a^l - ^\og\bi - bj\ - p loge 

\ i,j i<j i<j / 

(73) 
From this entropy function, we find infinite divisibility for the matrices (e '^'^(^i^j)) and hence 
conditional positivity of —SiyiVj). This follows from reflection positivity, since we can write 
for any A > 

(74) 
in terms of vertex (exponential) operators constructed with a free massless scalar field 0(x) 
|38j . The right hand side gives the left hand one since we have 

/Q|g*/'ia::/(a;)<^(a;)|g\ _ g^ / dxdj//(a;) log |z-y|/(j/) ^ /yg\ 

For general free fields the Renyi entropies are given in terms of exponentials of operators 
which are quadratic in the free fields. In this particular case, due to the bosonization 0741) . 
they can also be written in terms of exponentials of operators linear in a free scalar field. The 
infinite divisibility property coincides with the one of these exponential operators as discussed 
in Section 2 Q 

In this example, the Renyi entropies Sn are proportional to the entropy, Sn = ll^^)'^' ^^^ 
hence —Sn is conditionally positive. It is interesting to note that in the massive case the Renyi 
entropies can also be mapped to correlators of exponentials of a scalar field [3B]. However, 
this belongs to the interacting Sine-Gordon theory. Therefore, the infinite divisibility of the 
two point function discussed in the previous subsection corresponds to the one of correlators of 
these exponential operators in the Sine Gordon model. 

4.4 Two intervals in a CFT 



Formula fj73|) does not hold for general conformal field theories. Even if no exact result for the 
entropy is known for more than on interval in interacting models there are several interesting 
exact results for the Renyi entropies of integer n and two intervals for compactified scalars or 
the Ising model [HI |121 US, HI]. Most importantly, the exact results for the Renyi entropies 
allow us to show that in general the trp'^ are not infinite divisible in QFT [8] . An approximation 
expansion for the entropy in these models are found in [l5]. However, these approximations 
seem not to be suitable to test positivity. 

Conformal invariance implies the entropies of two intervals can be written as 

where Fn{x) = F„(l — x), F{0) = 1, is a function of the cross ratio x = (^Za)(b^b) > ^ is a 
constant, and c is the Virasoro central charge. The functions Fn{x) for the known examples 
are given by expressions involving theta functions and their inverses. A simple example is 5*2 
for the critical Ising model which is given in terms of algebraic functions |12] 



^^<^' ^ 7=2 



1/2 



^1 + yi)(l + y/l~Xj ^ , ^1/4 _ 1/4 , ^1/4 



+ x^/^ + (x(l-x))^/^ + (l-a;) 



1/2 



(77) 



■^The hypothetical theories with extensive mutual information considered in j40) also have infinite divisible 
-AS given by correlation functions of exponentials of free fields. 
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We have tested conditional positivity inequalities numerically and find that this Renyi entropy 
is not conditional positive. However, as expected, the non-linear inequalities ([6]) hold, what is 
consistent with trpy being given by a vacuum expectation value of a twisting operator. 

We do not know if the infinite divisibility is recovered in the so far unknown limit n — )■ 1 for 
the entanglement entropy. Since the Renyi entropies do not satisfy the inequalities, this case 
would be a very important test. 

4.5 Long distance limit for two regions 

In general, for free theories and two sets A and B which are far apart we have the mutual 
information is proportional to the squared of the Dirac or scalar field correlators [32] • Thus, 
this gives conditional positivity when only the distance d between A and B are changed but not 
the shapes of the sets. For the case of two intervals in a massive Dirac field in two dimensions 
we know also the dependence on the shape for large separating distances 



1 1 

I(A, B) ~ -m^ia+h^ + a^b^K^(md) + -m'^(a+b+ + a_b)Kf(md) , (78) 

6 6 

where a+, 6+, a_ and 6_ are the projections on the null coordinates of the intervals A and B, 
which are not assumed to lie in a single spatial line but can also be boosted to each other. The 
positivity of this quantity is proved by noting that the relevant matrix 

/ Kfimd) Klimd) \ 

\ Kl{md) Kl{md) ) ^ > 

is the Hadamard square of the correlation function of a massive Dirac field z7*^7^(\E'(x)\E'(y)^). 
This last correlator satisfies the Minkowski space analog of reflection positivity [32] ■ 

4.6 Angular sectors in 2+1 dimension 

In (2 + 1) dimensions, we study the infinite divisibility property in the entropy of planar angular 
sectors of a fixed side length L, for free scalar and Dirac fields. We place the angular sectors 
with common vertex in a plane xq^xi. We have for the entropy 

S{e) = {2 + 9)- + s{9) log(e) + finite . (80) 

The only divergent term in the entropy which does not cancel in the conditional positive 
inequalities is the logarithmic divergent contribution of the vertex angles. This has to satisfy 
the inequalities independently of the other finite contributions. The problem has then only one 
parameter. 

The checks require to solve with enough precision a coupled set of non linear differential 
equations and integrate the results on one parameter to obtain the logarithmic coefficient of 
the entropy [1^. Following ^\, this is done expanding the involved functions and differential 
equations in Taylor series around the value x = tt up to order (x — tt)^^. (for details see |47] ) . 
Again, as in the one dimensional case, up to what we have checked, we found the entropy 
satisfies the conditional positivity inequalities (158|) . as functions of the angle 6 of the angular 
sector, giving support to infinite divisibility for the exponential of the entropy. The figure (|5]) 
shows the functions /„ up to n = 5 (6 x 6 determinants), obtained from s{6) corresponding to 
a massless scalar. For the fermion field we find completely analogous results. 
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Figure 5: The logarithm of the functions /i, /2, /s, f^, /s (from top to bottom) obtained 
from the angular entropy function given by the coefficient of the logarithmic term s{6) of a 
plane angular sector for a massless scalar field. We evaluate these functions of the angle 6 as 
an expansion around 6* = vr up to order (vr — 6')^^ for s{6). These functions are positive in a 
neighborhood of 6' = tt. For small values of 6 the Taylor expansion for /„ does not give a good 
approximation, and this is the reason the last three functions go down steeply at some 6. 



4.7 Dimensional reduction for free fields 

In higher dimensions, for free fields, some universal terms in the entanglement entropy can be 
obtained via dimensional reduction from results calculated in lower dimensions [32] • Let us 
consider sets in k + d spatial dimensions of the form V = B x X, where B is a box on the first 
k coordinates xi,...Xk, of sides lengths Ri, i = l,...,k, and X is a set in d dimensions. The 
entropy of V in the limit of large Ri is extensive in the sides Ri. Thus, we can compactify the 
directions Xi, with i = 1, ...,k,hj imposing periodic boundary conditions Xi = Xi + Ri, without 
changing the result of the leading extensive term. In the limit of large Ri, the entanglement 
entropy S{V) is written in terms of S{X,m), the lower dimensional free entropy function for 
mass m. We have [35] 



S{V) 



kA 



2H^I'^T{k/2 



dpp^-'S{X, 



m^ 



pz 



where the transversal area A = Yli=iRi- The universal terms in S{V) will then come from 



the ones of X after this integration over the mass. From flSTl) . it is evident that if S{X,m) is 
conditionally positive, this is inherited by S{V) for these special type of regions. For example, 
the conditional positivity of the entropy of an interval in one spatial dimension is inherited by 
the entropy of an infinite strip in two dimensions. 



5 Some checks of the geometric inequalities for minimal 
surfaces 

In this section we check the inequalities which arise from conditional positivity of —S{V), using 
the holographic prescription (Ryu-Takayanagi ansatz) for the entropy in terms of minimal areas 
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in the bulk AdS space. We find counterexamples for the most complex geometries available, 
which depend on two parameters. Remarkable, the one parameter functions analyzed are 
conditionally positive. 

5.1 Long strips 

First consider the case where ^ is a rectangular strip of width R and length L, in the limit 
of large L, in three Euclidean dimensions. Because the result has to be extensive in L, and 
because of conformal invariance, we have 

S{V) = co--c^^, (82) 

for some constants Cq > 0, Ci, and a distance cutoff e. 

In order to obtain the positivity inequalities for single component regions of this kind, we 
position different regions in the x° > space, with the long sides lying parallel to the x^ 
direction. Further, we require one of the vertices to lie at the origin x = (see figure [6]) and all 
strips to have the same length L and orientation. Hence, all regions have a common boundary 
line. Because of the cancellation of the boundary lines passing on the point x = belonging to 
the regions and their reflected counterparts, the area is a two point function depending on the 
coordinates x = {x^, x^) of the free boundary line of each region. We now prove the conditional 
positivity of minus the function fl82|) . where R = \x — y\. Notice the conditional positivity 
condition does not take into account the cutoff dependent term in 0821) . 

We can expect the positivity of R~^, and hence the conditional positivity of —S{V), because 
R~^ is a field correlator in two dimensions. This requires ci > only. The same happens for 
any power R"" with u > 0, since R'" are field correlators for any u > 0: these are two point 
correlators of fields with different conformal weight in two dimensional CFT. It is instructive 
to go through the explicit proof in this two dimensional case. We are going to show that the 
correlators R~^ are positive for all i/ > 0. This is equivalent to — log(-R) being conditionally 
positive. This writes 



d'x £y fix)f{yr log((x° + y^ + {x' - y'f) > , (83) 

conditioned to J d'^x f{x) = 0. Writing the correlators as a Laplace-Fourier transform 

log((x° + yy + {x' - y^f)) = I dko dh e'^'^^'-y'^-^'^^^'+y'^g{h. h) , (84) 

and replacing this in (18 3 p we have 

dkodh\fik)\^giko,h)>0, (85) 



with 

/>) = / d'x 6-'-^°+'"^^' fix) , /(O) = 0. (86) 



This last Laplace transform makes sense since the condition for reflection positivity is f{x) = 
for x^ < 0. Then, all the infinite series of inequalities are summarized as g{ko, ki) > for 
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Figure 6: Three rectangular regions Vi, V2 and V3, and their time reflected images, with common 
side at the origin. Only the short sides of the rectangles are shown. The loop V1V2 is equivalent 
to a rectangular loop with side R, in the limit of large longitudinal size L ^ R. 

k ^ 0. The explicit evaluation of the Laplace anti-transform of the logarithm confirms this 

expectation 

., , s S{ko- \h\) 
g{ko, fci) = — . (87) 

ml 

For higher dimensions there is a minimal power u such that R~'^ is a positive correlation 
function. This unitarity bound for two point functions in CFT is u > {d — 2). The proof can 
be done using the Laplace transform as above. Notice that for dimensions d > 2 this means 
there is a lowest power for which the two-point correlator can be positive and hence there are 
no infinite divisible two-point correlation functions. However, this is not the case for extended 
objects such as surface operators of surface dimension q. Then in the conformal case we expect 
correlators proportional to e^'^^'', where L is the large size of the branes in the q directions, and 
R is the separation between two branes. Thus R lives in d — q dimensions. The unitarity bound 
implies this correlator is infinite divisible {R~'^ is conditionally positive) only ii q > {d — 2)/2. 
For the case of the entropy q = d — 2 and we always have infinite divisibility of strips in any 
dimension. This is also the case of Wilson loops for free electromagnetic field in d = 4, which 
has q = I- 



5.2 Angular sectors 

Now consider loops which bound a plane angular sector of angle 6. The area as a function of 
the angle is [M] 

A{C) = ci + C2(2 + 9)-- g{e) log(L/e) . (88) 

e 

The function g{6) is given in parametric form in terms of a variable x G (0, 00) as 

dz 



9 



2xVr 



x^ 



(^2 + a;2)^(22 + x2 + 1)(22 + 2x2 + 1) ' 



V Vz^ + 2x2 + 1 



(89) 
(90) 
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Figure 7: The logarithm of the functions fi{6), with i = 1 to i = 5 (the curves with higher 
slope correspond to larger i). These functions are all positive (the logarithms are real). They 
are symmetric around 9 = tt, fi{6) = fi{2n — 6). 



This gives g{9) for 9 between ^ = (x — )■ oo) and ^ = tt (a: — t- 0). We can extend it to ^ G [0, 27r] 
by using ^(27r - ^) = g{9). 

We consider angular sectors on the plane (x'^,x^), of the same size L with common vertex 
at x = 0, and common side along the x^ axes. Then, the first two terms in (I88p are not relevant 
for the inequalities. These write 



d9^ / d92 u{9^)g{9i + ^2)^(^2)* > , 



(91) 



for arbitrary u{9) with support in G [0,-7r], and J d9u{9) = 0. 

If we could write g{9) as a Laplace transform (perhaps for a distributional or discrete 
function g{p)), 

g(9)= [dpe~P'g{p), (92) 



we could simplify the inequalities as g{p) > using the same trick as in the previous section. 
Unfortunately at present we cannot understand whether an expression like f l92|) is valid. Thus, 
we can not offer a complete analytical check of these inequalities here. 

However, we have checked numerically the functions fn{9) of the derivatives oi g{9) which are 
predicted to be positive by the inequalities flSSl) . Up to what we have checked, these functions 
are indeed positive (see figure [7]). We have also checked the positivity of the determinants using 
up to 6 random angles entropy matrices. 



5.3 Coplanar circles 

We consider two coplanar circles of radius -Ri and -R2; at a distance h = hi + h2, where hi and 
/i2 are the distances of the circles to the plane x^ = 0. This case contains objects determined 
by two parameters, {h, R). The area for two parallel coaxial circles (rather than coplanar) was 
calculated in [19] and [50]. This case can be transformed conformally to the one of the two 
coplanar circles [SI]. Taking out a divergent piece proportional to the circles perimeter, the 
relevant part of the area is conformally invariant. It is then given as a function of the cross 
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Figure 8: The logarithm of the functions /i(-z) for i = 1...8 (from bottom to top) corresponding 
to coplanar circles of equal radius (here Ri = R2 = 1). They are all positive in the range 
z e (0,2:0), corresponding to connected minimal surfaces. 

ratio z = , ~^"^ , "^'^ of four points. We choose xi and xa to be the points in the two circles 
which are at a greatest distance to each other, and X2 and X3 the ones at the minimal distance. 
Then the cross ratio is 2 = h- — |^-^ — - E (0, 00). We have a parametric expression for the 

cirGcL 
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asin2(t) + a/1 + asin2(t)) 
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(93) 



(94) 



where E and K are the complete elliptic integrals of the first kind. The parameter x is related 
to the geometry of the configuration by the following equation involving the cross ratio z 



1 / \ /"^ 

- log fl + 2^ + 2^fz^T~z\ = X / 

V2x [k (- 



, / V43:2+l-l 



sin2(0) 



a/cos2((/)) — x2sin^((/)) 



2x2 



n 



l+Vl+4x2 



2x2 



2x2 



(95) 



1 + 4x2-1 

where n(x|y) is the complete elliptic integral of the third kind. 

We are interested in the solution which describes a surface connecting the two loops. For 
each z < zq = 0.27288553... there are two solutions of equation (jHSD for x, but only one 
gives the minimal area. This corresponds to the solution with the greater value of x (having 
X > 0.58110028...). The other solution corresponds to loops having the same orientation, which 
is not the case appearing in the reflection positivity inequalities (and it does not satisfy these 
inequalities). 

Considering the case of circles of equal radius R\ = R2 we have that the entropy is a one 
dimensional function and we evaluate the corresponding functions /„. These are shown in flgure 
dH]). This one parameter case seems to lead to a conditional positive function. 
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Figure 9: The logarithm of the absolute value of the eigenvalues of the matrix of infinitesimal 
inequalities of order 5 as a function of z (there are global factors depending on R which we fix 
by setting Rx = R2 = V). One of the eigenvalues becomes negative for z > 0.165. 




Figure 10: The logarithm of the absolute value of the eigenvalues of the matrix of infinitesimal 
inequalities of order 9 as a function of z (there are global factors depending on R which we fix 
by setting Ri = R2 = I). There are negative eigenvalues for z > 0.09. 
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Then, we consider the full two parameter case, and evaluate the eigenvalues of the matrices 
in the derivatives discussed in Section 4 with respect to both parameters, radius R and distance 
h to x^ = plane. The first non trivial matrix (1^2]) has 5 eigenvalues which we evaluate at 
the point of equal radius Ri = R2. We plot these eigenvalues as a function of z in figure 1^. 
We see this has a negative eigenvalue for z > 0.165. Figure flTOj) shows the eigenvalues of the 
matrices at the next order which consists of 9 x 9 matrices. Negative eigenvalues appear now for 
z > 0.09. According to the general theory, since the functions are definitely not conditionally 
positive in this case, the matrices of higher rank should be able to detect negative eigenvalues 
for z as small as we want. 

The same analysis can be carried out for the case of two parallel concentric circular Wilson 



loops. The function of the cross ratio is the same, but this is given as z = ^"'" ar r ^ — 



4/?l/?2 



G 



(0, 00) in terms of the radius -Ri and R2 of the circles, and their distances hi and /i2 to the 
x° = plane (which in this case is assumed parallel to the circles). The results are completely 
analogous to the ones above. In fact it can be seen that a conformal transformation commuting 
with the time reflection operation maps the positivity relations in one case and the other. 
Hence, we have infinite divisibility for the one parameter case of circles with equal radius, but 
not for the full two parameter function. 

5.4 Two spheres in four dimensions 

The case of two spheres can be treated similarly. Again the relevant part of the minimal area 
is conformally invariant and depends on a cross ratio z = h- — j^-^ — - G (0, 00). The formula 
for the area can be obtained from the result for two concentric shells [22] by a conformal 
transformation. We have for the mutual information of the two spheres A and B, 

sin2(t0) 



, pn«2|'/„^ rsm^m 

log(l + 2z + 2v/^T^) = ^^ / dy^= 

sm'^(to) Jo J(l 



I(A,B) = lim 

€->0 



v/(l-2/)((l-2/)2-2/3cos4(to)/sin«(to)' 



*o , cos^rt) r*" , cos^(t) 

dt ^ ^ - / dt ^ ' 



.3m . Lr.^4:(+\ _ ^^.AU \ sm^(*) J^ ^\^'i(+\ . nr.^4:(+\ _ o^o4/'+ N sin6(t) 



„, sm\t)^cos\t) - cos4(to) J^ -^1% sin^(t) Y^cos4(t) - cos4(t,,^.^,(^^^_, 

The connected minimal surface exists for z < 0.0965. 

We obtain results very similar to the ones of the previous subsection. The one parameter 
function for equal radius -Ri = R2 gives place to (apparently) a conditional positive function. 
When the spheres are also allowed to vary in size the conditional positivity is lost, and we find 
negative eigenvalues for the matrices of multidimensional infinitesimal inequalities. We do not 
present the relevant figures here since they are in all respect similar to the ones of the last 
subsection for circular regions. 

5.5 Multicomponent strips and the topology of the surfaces 

Now we consider regions formed by an arbitrary set of disjoint parallel strips. The analysis 
of infinite divisibility will display a curious phenomena, related to the Gross-Ooguri phase 
transition [52]. 
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We take all long edges on the regions parallel to the x^ direction, with size L. Let us 
call ai...an to the lines on the boundary with one orientation and bi...bn to the lines with the 
opposite orientation, where n is the number of connected components. In the large L limit the 
relevant piece of the minimal surface will be some choice of all the possible surfaces formed by 
the union of the minimal surfaces connecting just two lines, one of type a and one of type b. 
This can then be written as (ai, 6o-(i), ^2, &cr(2), •••, On, ba{n)), with <j{i), the permutation of l...n, 
producing the minimal area. This minimal area is just the sum A{ai, fe(j(i)) + ... + A{an, &o-(n))- 

Let us then think in a correlator involving m of these multicomponent Wilson loops Vi, 
i = l...m. We are interested in the matrix A{ViVj). Let Vi have 2nj boundary lines. Then 
the minimal area of ViVj contains a group of 2qij < 2ni surfaces which cross the x" = plane 
and are connected to a group of 2qij lines in Vi and 2qij lines in Vj. There is also a group of 
2ni — 2qij lines in Vi which are connected among themselves, decoupled from Vj. Then, suppose 
that in any of the correlators involving Vi in the matrix, the same group of lines of Vi are 
connected among themselves. The contribution of these lines is necessarily the same in all the 
correlators involving Vi in the correlator matrix in question. Let us call this contribution C,i- 
Now suppose the same happens for all the regions Vk, that is, all the minimal surfaces in the 
correlator matrix elements VkVj which contain Vk always has the same subset of the lines of Vk 
which are connected among themselves. Then we write 

A{V^^,)=A{V^^,)+^, + ^,, (97) 

where A{VVj) contains the contribution from the surfaces which cross the x° = plane. Eq. 
fpUj) then tells that —A{VVj) is conditionally positive if and only if —A{ViVj) is conditionally 
positive. In this situation, the number of surfaces which cross x° must be the same in all 
correlator in the matrix, and hence g^ = g is independent of i,j. In this case the problem is 
simplified, we can forget about the disconnected pieces and think all lines are connected across 
x° = by the minimal surface, and all loops have the same number of lines 2q. 

Once we have disposed off the disconnected pieces, we prove infinite divisibility under the 
following conditions: For all i,j we require there is a smooth transformation Vij{t), for t G [0, 1], 
such that Vij{0) = Vi, Vij{l) = Vj. This requires in particular that all the regions have the 
same topology. Further, we ask these transformations to form a group under concatenation 
of mappings, Vij o Vjk = Vik. And finally, but most importantly, we require that the sequence 
of minimal surfaces corresponding to Vii'{t)Vjji{t), going from the minimal surface of VVj to 
the one of Vi'Vj', to be smooth. That is, the trajectory of minimal surfaces does not have to 
cross any phase transitions in the middle. In particular, the minimal surfaces on the correlator 
matrix must belong the same homotopy class of minimal surfaces. 

It is immediate that this condition implies for the kind of loops we are considering here that 
the minimal surfaces connect only homologous lines. More clearly, let us name the lines in Vi 
as a^ , ..., ttq and b^ , ..., bq . The conditions above imply we can choose the names of the lines 
such that in the minimal surface of ViVj the line a^ is connected with a)^ , and b^ with b^ , for 
all i,j, k. With this naming of lines, all the minimal surfaces of the correlator matrix connect 
the k^^ lines of type a, and b among themselves. Then we can write 



- A{VV,) = - J2{A{a'i\ a^'^) + A{l^,V^^)) . {91 



k=l 

This is a sum of conditionally positive matrices and then conditional positive. 
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Conditional positivity cannot be proved in the case where the surfaces involved in the 
correlator matrix are separated by phase transitions, and some numerical examples show that 
generically it does not hold in these cases. 

Then, we can think that for more general cases than the one studied here the conditional 
positivity, except for the disconnected components, requires the minimal surfaces in the corre- 
lator matrix to lie in the same sector in the sense specified above. Notice that this does not 
exclude any topology for the minimal surfaces involved in the correlator. Generally, the con- 
ditions can always be achieved if the loops participating in the correlation matrix lie in some 
sufficiently small neighborhood of a given loop V. However, even in this cases without phase 
transitions we have already found some counterexamples of the inequalities in section 5.3. 

6 Final comments 

Motivated by the positivity relations satisfied by the integer n index Renyi entropies, and a 
natural conjecture about the path integral representation for the n — >■ 1 limit, we have checked 
whether the entanglement entropy gives place to conditionally positive correlation functions. 
We could not find counterexamples of the inequalities ([9j) among the few known exact QFT 
results. This is not so for the minimal areas of the holographic entanglement entropies where 
we found counterexamples for the case of two spheres or circles of different radius. In this 
case, the two parameter infinitesimal inequalities fail. The relation between analyticity and 
positivity then implies that the inequalities have to fail in the whole domain of holomorphy of 
the two variable correlator. We have also found phase transitions are a source of violations of 
the inequalities for multicomponent regions. 

A question which remains open is the status of the conditional positivity inequalities for 
the entropy in QFT. In the light of the results for minimal surfaces, the evidence in favor 
of the inequalities in QFT is quite weak. On one hand, most of the discussed cases refer to 
entanglement entropies for free theories, which have there own special properties and relations 
to the Renyi entropies. On the other, for the simple geometries analyzed, the minimal areas also 
obey the inequalities. The entanglement entropy for two spheres is not known exactly for a QFT 
yet. Of course, a positive answer for the infinite divisible property of the entanglement entropy 
in QFT would mean the holographic entropy ansatz has to be modified for some geometries. 
Perhaps the full solution of the entanglement entropy for a simple geometry such as the case 
of two intervals in general CFT would be enough to settle this question. 

Another interesting question is the direct mechanism (either in the context of QFT or a 
purely geometric one) behind the validity of the inequalities in the one-parameter functions 
discussed in this paper. Indeed, we have found evidence indicating that several non trivial 
functions are conditional positive, obeying an infinite series of inequalities on the function 
derivatives. This includes minus the entropies of one interval and one angular sector for free 
scalar and fermions, and the minimal areas corresponding to angular sectors, and pairs of circles 
and spheres of equal radius. In terms of the minimal areas the conditional positivity property 
is rather peculiar. For example it is easy to check it does not hold for ordinary geodesies in 
AdS, but only for the limit of geodesies starting and ending at the boundary. 

One possible explanation could reside in the interplay between the positivity recurrence 
relations and some unknown differential recurrence relations for these functions. On the minimal 
surfaces these examples of infinite divisibility may be hinting to the validity for the small 
angle path integral representation of the entropy in these geometries. It would be particularly 
interesting to understand and generalize this enlargement of cases beyond the single sphere and 
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the plane. 

It has been discussed in the hterature that the small angle path integral representation is a 
weakness in Fursaev's proof of Ryu-Takayanagi proposal for holographic entanglement entropy 
[TT| m] . We have further found it seems to involve contradiction for certain geometries. On 
the QFT side the assumptions on which the proof is based lead to infinite divisibility, but 
the holographic result is in general not infinite divisible. However, for some geometries this 
contradiction is absent. 
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